CRITICAL WINDOW FOR THE VACANT SET LEFT BY RANDOM 
WALK ON RANDOM REGULAR GRAPHS 
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Abstract. We consider the simple random walk on a random d-regular graph with n 
vertices, and investigate percolative properties of the set of vertices not visit ed by the walk 
until time un, where u > is a fixed positive parameter. It was shown in jCTWll] that 
this so-called vacant set exhibits a phase transition at u = it*: there is a giant component 
if u < u+ and only small components when u > u*. In this paper we show the existence of 
a critical window of size n — 1 / 3 around u*. In this window the size of the largest cluster is 
of order n 2 / 3 . 

o 

03 1. Introduction 

The study of percolative properties of the vacant set left by a random walk on finite 
i graphs was initiated by Benjamini and Sznitman |BS08j for the case of random walk on a 

, , high-dimensional discrete torus (Z/iVZ) '. In [BS08] it is proved that if the random walk 

runs up to time uN d , where u is a small constant, the vacant set has a giant component 
with volume of order N d , asymptotically as N grows. On the other hand, if instead we 

^ consider a large constant u, all components of the vacant set have a volume of order at most 
log A N (for some A > 0) as was proved in [TW10J. This shows the existence of two distinct 
phases for the connectivity of the vacant set on the torus as u varies. However, the above 
mentioned works leave several open questions, such as whether the transition between these 

t-h two phases happens abruptly at a given critical threshold and, if this is the case, how does 
the vacant set behave at the critical point? 

Such problems are much better understood when instead of the torus one considers ran- 
dom ci-regular graphs or <i-regular large-girth expanders on n vertices. In this case, when 
random w alks runs up to time un, it is known that the vacant set exhibits a sharp phase 
transition (CTWllj : the size of the largest connected component of the vacant set drops 
abruptly from order n to order \ogn at a computable critical value «*. In this paper we 
explore more closely this phase transition, in particular we prove that the size of the largest 
component of the vacant set exhibits a double-jump, similar to that observed in Erdos-Renyi 

^ random graphs. 

Let us now give a precise definition of the model. Let Q n ^ be the set of all non-oriented 
o?-regular simple graphs with n vertices (here and later we tacitly assume that nd is even). 
Let F nt d be the uniform probability distribution on Q n ,d- For any graph G = (V,S) let P G 
denote the canonical law on the Skorokhod space D([Q, oo), V) of a continuous-time simple 
random walk on G started from the uniform distribution. We use (X t )t>o to denote the 
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canonical coordinate process. For a fixed parameter u > 0, we define the vacant set as the 
set of all vertices not visited by the random walk before the time u\V\, 

(1.1) V u = V G = V\ {X t : < t < u\V\}. 

We use C max to denote the maximal connected component of the subgraph of G induced 
by V". The vacant set and in particular its maximal connected component are the main 
objects of invest igation in this paper. 

As proved in CTWllj . the phase transition in the connectivity of the vacant set occurs 
at the value u* given by 

d(d-i)\n(d-i) 

(L2) u * = — {d^w — ' 

and can be described as follows: Let G n be a graph distributed according to P njC j. Then 
with P ra! d-probability tending to one as n — > oo: 

Super-critical phase: For any u <u± and o > there exist p and c depending on u, 
a, and d, such that 

(1.3) P G "[\Cl.J>pn]>l-cn-°. 

Sub-critical phase: For any u > u± and a > there exist K and c depending on u, 
cr, and d, such that 

(1.4) P G ~[\Cl.J<K\ogn]>l-cn-°. 

In this paper we study the behaviour of the vacant set in the vicinity of the critical point. 
The main results of this paper are the following two theorems. In their statement we use 
P n ,d to denote the annealed measure 



(1.5) P M (-) = J P G (-)P M (dG). 

We say that an event A occurs P ni( rasymptotically almost surely (or simply P n d -a.a.s.) if 
lim^oo P n ,d( A ) = !• 

Theorem 1.1 (Critical window). Let (u n ) n >i be a sequence satisfying 

(1.6) \n l ^{u n — u±)\ < A < oo for all n large enough. 

Then for every e > there exists A = A(e, d, A) such that for all n large enough 

(1.7) VnM' 1 ^ < ICLl < An 2 / 3 ] >l-s. 

If u n is not in the critical window, then the maximal connected component behaves 
differently: 

Theorem 1.2. (a) When (u n ) n >i satisfies 

(1.8) u+ — u n > U, ana u n := n 1 (u* — u n ) > oo, 

then for v n = 2n 2 ^ 3 u n ^■-^ 2 e - u *i d - 2 )/ ( . d - 1 ) an d for every e > 



\PU n I 

I '-'max I 



< e P n , d -a.a.s. 



1.9) 

(b) When {u n ) n >i satisfies 
T.1UJ u* — u n > U, ana u n := n 1 (u* — u n ) > —oo, 



then for every e > there exists B = B(e) > 0, such that for all n large enough 



;i.ii) 



nA 



max 



< Bn 2 '* 



-1/2 



]>i 



The above theorems confirm that the maximal connected component of the vacant set 
behaves similarly as the largest connected cluster of the Bernoulli percolation on random 
regular graphs, see [ABS04|, INP10| IPit08j . Remark that the upper bound in part (b) of 
Theorem L2 seems to be non-optimal. The result is however sufficient to confirm that the 
width of the window is n -1 / 3 . To improve such a statement, it is necessary to obtain better 
estimate in Theorem |3.2|(ii) that we quote below. 



The methods of this paper are largely inspired by the recent article by Cooper and Frieze 



|CF10j . where the authors develop a new technique to prove (1.3) and (1.4). This technique 
is very specific to deal with random regular graphs, in contrast with the results of |CTWll| 
which hold for a more general class of graphs, including e.g. large girth expanders. In 
the present article we extend the methods in |CF10| to the critical case. We believe that 
obtaining such an extension from the techniques in CTW11 should be rather difficult. 

The crucial observation of |CF10] . allowing for a very elegant proof of (1.3), (1.4) for 
random regular graphs is the following. Under the annealed measure ( 1.5[ ), given the infor- 
mation about the graph discovered by the random walk up to time u\V\, the subgraph of 
G induced by the vacant set V u is distributed as a random graph uniformly chosen within 
the set of graphs with a given (random) degree sequence, see Proposition 3.1 below. 



The paper |CF10j further uses the fact that the behaviour of the uniform random graphs 
Gd with a given degree sequence d : V — > N is sufficiently well known. More precisely, as 



follows from |MR95] . there exists a single parameter Q = Q(d) (see (3.6) below) such that 
Q < implies that Gd is typically sub-critical (i.e. has only small components), and Q > 
implies that Gd is supercritical (i.e. has a giant component). Moreover, the recent paper 
[HMIO] establishes the existence of an intermediate regime (the so-called critical window) 
when Q(d) converges to zero at a certain rate as n tends to infinity, see also |JL09j . The 
results mentioned in this paragraph that will be useful in this paper are summarised in 
Theorem 13.21 below. 

The principal contribution of this paper is thus to obtain sufficiently sharp estimates on 
the random degree sequence of the vacant set V", and consequently on the value of the 
parameter Q, see Theorems 5.1, 5.2 and (6.3) below. Weaker estimates of this type were 



shown in |CF10j . which combined with |MR95] . allowed them to deduce (1.3), ( |1.4 ) 



We should remark that [C FlOj contains also a statement on the critical behaviour. More 
precisely, Theorem 2(iii) of |CF10] states that for some u n = + o(l)) (which might 
be random), the size of C^™ x is n 2//3+ °^\ P n ^-a.a.s. Our results considerably improve this 
statement. 

Note also that much more is known about the random graphs with a given degree sequence, 
see for instance the results in |FR09] and |JL09j . Often the hypothesis of these results can be 
shown to hold true for V u , using Theorems 5^ and 5J2 If this is the case, their conclusions 
will also apply to V u (P n) d-a.a.s) 
the vacant set. 



providing us with more information on the geometry of 
As an example of such application, we obtain the following improvement on the statement 



(1.3) about the super-critical behaviour of the vacant set. 



Theorem 1.3. Let u < it*. Then there is p = p(u, d) G (0, 1) such that for every e > 
(1.12) n-^C^Jeip-^p + e) P n4 -<,.a.s. 



e(p-e,p + e) 
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In the above statement, the value of p can be explicitly calculated, see (6.9 ) below. Remark 
also that to obtain the above theorem, our precise estimates on Q are not necessary, in fact 
the precision obtained in |CF10] would have been sufficient. 

The former, establishes an estimate 



Finally, let us briefly describe Theorems 5.1 and 5.2 



on the expected degree distribution of V u , by approximating the probability that a random 
walk visits a neighbourhood of a given vertex x G V before time un. For this, we make use 
of the well known relation between random walks on graphs and discrete potential theory, as 
well as the pairing construction introduced by Bollobas, which we detail in Section [2] Then 
in Theorem 5.2 we prove that with high probability the degree sequence of V" concentrates 
around its expectation. This is done using a standard concentration inequality, together 
with the fast mixing properties of the random walk on a random regular graph. 

This paper is organised as follows. In Section [2] we introduce some of the notation needed 
in the paper and the pairing construction of random regular graphs. In Section [3j we recall 
the results of |GF10[ IHM10[ IJL09] needed later. Section [1] contains precise estimates on 
the behaviour of the simple random walk on random regular graphs. In Section [5| we give 
the estimates on the degree sequence of the vacant set. Theorems 1.1-1.3 are proved in 
Section [6j The Appendix summarises some general facts concerning random walks on finite 
graphs. 



2. Notation and definitions 

2.1. Basic notation. We now introduce some basic notation. Throughout the text c or 
d denote strictly positive constants only depending on d, with value changing from place 
to place. Dependence of constants on additional parameters appears in the notation. For 
instance c u denotes a positive constant depending on u and possibly on d. We write N = 
{0, 1, . . . } for the set of natural numbers, and [d] for the set {1, . . . , d}. For a set A we 
denote by \A\ its cardinality. For any sequence of probability measures P n and events A n 
we say that A n holds P n -a.a.s. (asymptotically almost surely), when lim^oo P n L4 n ] = 1. 

In this paper, the term graph stands for a finite simple graph, that is a graph without 
loops or multiple edges. Sometimes we intentionally allow the graph to have loops and/or 
multiple edges and in this case we use the term multigraph. For arbitrary (multi)graph 
G = (V, £), we use dist(-,-) to denote the usual graph distance and write B(x,r) for the 
closed ball centred at x with radius r, that is B(x,r) = {y G V : dist(x,y) < r}. 

We use Gn,d (resp. J^i n ,d) to denote the set of all ci-regular graphs (resp. multigraphs) 
with vertex set V n — {1, ... , n}. Given a degree sequence d : V n — > N, we use Gd to denote 
the set of graphs for which every vertex x G V n has the degree d x = d(x). Similarly, Aid 
stands for the set of such multigraphs; here loops are counted twice when considering the 
degree. P n ^ and denote the uniform distributions on G n ,d and Gd respectively. 

2.2. Pairing construction. In order to study properties of random regular graphs, Bol- 
lobas (see e.g. jBolOlj ) introduced the so-called pairing construction, which allows to gen- 
erate such graphs starting from a random pairing of a set with dn elements. The same 
construction can be used to generate a random graph chosen uniformly at random from 
Gd- Since this pairing construction will be important in what follows, we give here a short 
overview of it. 

From now on, whenever we consider a sequence d : V n — > N, we suppose that ^2 xeVn d x is 
even. Given such a sequence, we associate to every vertex x G V n , d x half-edges. The set of 
half-edges is denoted by Hd = {(x,i) : x G V n ,i G [d x ]}. We write H n ^ for the case d x = d 
for all x G V n . Every perfect matching M of Hd (i.e. partitioning of Hd into \H&\/2 disjoint 
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pairs) corresponds to a multigraph Gm = (Ki,£m) £ M-d with 

(2.1) E M = {{x,y} : {{x,i), {y,j)} G M for some i G [d x ], j G [d y ]}. 

We say that the matching M is simple, if the corresponding multigraph Gm is simple, that 
is Gm is a graph. With a slight abuse of notation, we write Pd for the uniform distribution 
on the set of all perfect matchings of Hd, and also for the induced distribution on the set 
of multigraphs Aid- It is well known (see e.g. [BolQl] or |McD98j ) that a Pd distributed 
multigraph G conditioned on being simple has distribution Pd, that is 

(2.2) f d [ Ge .\ Ge g d ] =¥d [ Ge .] ) 

and that, for d constant, there is c > such that for all n large enough 

(2.3) c < F n4 [G e g n4 ] < 1 - c. 

These two claims allow to deduce P n)( j-a.a.s. statements directly from P nj d-a.a.s. statements. 

The main advantage of dealing with matchings is that they can be constructed sequen- 
tially: To construct a uniformly distributed perfect matching of Hd one samples without 
replacements a sequence hi, ... , h\H d \ of elements of Hd in the following way. For i odd, hi 
can be chosen by an arbitrary rule (which might also depend on the previous (/i 3 -)i<i)> wm l e 
if i is even, hi must be chosen uniformly among the remaining half-edges. Then, for every 
1 < i < \Hd\/2 one matches h 2 i with hn-\- 

It is clear from the above construction that, conditionally on M' C M for a (partial) 
matching M' of Hd, M \ M' is distributed as a uniform perfect matching of Hd \ {{x,i) : 
(x,i) is matched in M'}. Since the law of the graph Gm does not depend on the labels %' 
of the half-edges, we obtain for all partial matchings M' of Hd 

(2.4) P d [G M \M' e-\MDM']= F d '[G M G ■}, 

where d' x is the number of half-edges incident to x in Hd that are not yet matched in M' , 
that is d' x = d x - \{{(yi,i),(y2,j)} G M' : yi = x,i e [d x ]}\, and G m \m> is the graph 
corresponding to a non-perfect matching M \ M', defined in the obvious way. 

2.3. Random walk notation. For an arbitrary multigraph G = (V,£), we use P x to 
denote the law of canonical continuous-time simple random walk on G started at x G V, 
that is of the Markov process with generator given by 

(2.5) Cf(x) = J2P*y(f(v) - /( x ))> for / : 1/ -)• R, i G y. 

y&V 

Here p xy = n xy /d x , n xy is the number of edges connecting x and y in G, and d x is the degree 
of x; the loops are counted twice in n xx and d x . 

We write P x ' e for the restriction of P x to D([0, £}, V) and P xy ,£ for the law of random walk 
bridge, that is for P x > e conditioned on Xt = y. We write E x , E x ' 1 , E xy ^ for the corresponding 
expectations. The canonical shifts on D([0, oo), V) are denoted by Qf The time of the n-th 
jump is denoted by r n , i.e. r = and for n > 1, r n = mi{t > : X t ^ X } o Q Tn _ x + r n _j. 
The process counting the number of jumps before time t is denoted by N t = sup{k : < t}. 
Note that, when G is simple, under P x , (N t ) t > is a Poisson process on M + with intensity 
one. We write X n , n G Z + , for the discrete skeleton of the process X t , that is X n = X Tn . 

Given A C V, we denote by if^ and i?^ the respective entrance and hitting time of A 

(2.6) H A = inf{t >0:X t e A}, and H A = H A o 9 Tl + n. 

We denote by n the stationary distribution for the simple random walk on G, which is 
uniform if G is <i-regular (even if G is not simple). P G stands for the law of the simple 
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random walk started at 7r and E G for the corresponding expectation. For all real valued 
functions /, g on V we define the Dirichlet form 

(2.7) V(f,g) = \J2 (/(*) ~ f(y))(9(x) - 9(v))**Pxv = -^Cf(x)g(x)ir a . 

x,y£V xeG 

The spectral gap of G is given by 

(2.8) X G = mm{V(f, f) : n(f 2 ) = 1, tt(/) = 0}. 
From |SC97j . p. 328, it follows that for <i-regular graphs, 

(2.9) sup \P x [X t = y\- n y \ < e~ XGt , for all t > 0. 

x,yEV 

It is also a well known fact (see e.g. |Fri08j ) that there exist a > such that 

(2.10) Xq > a, both P„ irf -a.a.s. and P^-a.a.s. 

3. Preliminaries 

3.1. Distribution of the vacant set. Recall the notation V u = Vq for the vacant set of 



the random walk on the graph G = (V,£) at level u, (1.1). For the purpose of this paper, 
it is suitable to define a closely related object, the vacant graph V" = (V,£ u ) where 

(3.1) £» = {{x,y}e£:x,yeV%}. 

It is important to notice that the vertex set of V" is V and not V u , in particular V" is not 
the graph induced by V u in G. Observe however that the maximal connected component 
of the vacant set C max (defined before in terms of the graph induced by V u in G) coincides 
with the maximal connected component of the vacant graph V u (except when V u is empty, 
but this difference can be ignored in our investigations). 

We use T> u : V — > N to denote the (random) degree sequence of V u , and write Q^ d for 
the distribution of this sequence under the annealed measure P n d, defined by P n d(") : = 

The following proposition from Cooper and Frieze |CF10j allows us to reduce questions 
on the properties of the vacant set V u of the random walk on random regular graphs to 
questions on random graphs with given degree sequences. 

Proposition 3.1 (Lemma 6 of |CF10j ). For every u > 0, the distribution of the vacant 
graph \ u under is given by where d is sampled according to d , that is 

(3.2) P M [V"e-] = J f d [G e -}Q u n , d (dd). 



Although a proof of Lemma 3.1 can be found in |CF10j . we provide a proof here for the 
sake of completeness. 

Proof. Let M be a P niC r distributed pairing of H n d and let X be a random walk on G = Gm- 
Define M t C M to be the set of all pairs of half-edges incident to a vertex X s with s <t, 

(3.3) M t ={{(x,i),(y,j)}eM:xe{X s :s<t},ie[d}}. 

It is easy to see that the edges of the vacant graph V correspond exactly to the pairs in 
M\M un , that is V" = Gm\m„„- In particular, V u (x) is the number of the half-edges incident 
to x not matched in M un . Denoting by T u the a-algebra generated by ((X s , M s ), s < un), 
the above implies that T> u is J-^-measurable. 
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It follows from (2.4) that, conditionally on T w the distribution of Gu\M un ov ^f depends 
on the sequence of half-edges that are not matched in M un , and is given by Pd u . More 
precisely, 

(3.4) P M [V" e ■ \T U ] = P n , d [G MWun e ■ \V U ) =P B .[Ge-] ) 
and thus 

(3.5) P M [V U G ■] = P M [P M [V W G • = P M [iV[G G •]] = J ¥ d [G G • iQJ^dd), 

where the last equality follows from the definition of d . This concludes the proof of 
Proposition 3.1| □ 

3.2. Behaviour of random graphs with a given degree sequence. We now summarise 
the results about the behaviour of random graphs with a given degree sequence which will 
be used in this paper. For a degree sequence d : V n — > N we define 



(3.6) 



T n d 2 
Q(d) - ^ x=1 x 



and set rij(d) to be the number of x < n with d x = i, 
(3.7) n i (d) = \{xeV n :d x = i}\. 

For any graph G we use C max (G) to denote the maximal connected component of G. 
The following theorem summarises the results of [MR95, IJL09t IHM10] needed later. 

Theorem 3.2. Let (d ri ) n > 1 , d n : V n — > N, be a sequence of degree sequences. We assume 
that the degrees are uniformly bounded (i.e. max{d™ : n > 1, x < n} < A < oo), and that 
ni{d n ) > (n for a ( > 0. 

(i) (critical window) If \Q(d n )\ < AriT 1 / 3 for all n > 1, then for every e > there exists 
A = A((, X,e, A) such that for all n large enough 

IVL4-V/ 3 < \C max (G)\ < An 2 / 3 ] > 1 - e. 

(ii) (below the window) If lim n _ s . 00 n l l 3 Q(d n ) = — oo and lim^oo Q(d n ) = 0, then for 
every e > exists B = B((,e, A) < oo such that for all n large enough 



|C max (G)| <By/n/\Q(d 



(3. 



(iii) (above the window) Let lim n _ 
tion, assume that 

n t (d n ) 



n 



1/3 Q(d n ) 



> l-e. 



-oo and lim r 



Q(dT 



0. In addi- 



lim 

n— >oo 



n 



Pi 



for all < i < A, 



for some probability distribution (pj)o<i<A on {0, ...,A} ; and set X = ^f =Q iPi, 
(3 = J2t=oK^ ~ — 2)pi, and t> n = 2nX 2 [3~ 1 Q(d n ). Then, for every e 

|C ma x(C)| 



- 1 



< e, 



Pd"-a.a.s., 



(iv) (super- critical regime) Let \im n ^ 00 Q(d n ) = > and assume that (3.8) holds. 
Let g be the generating function of(pi), g(x) = 2~2t=oPi xl - Then there exists a unique 
solution £ to g'(x) = Xx in (0, 1), and for p = 1 — g(£) and any e > 



|^max(C) | 



< e, 



n 



Pd"-a.a.s. 



Proof. Parts (i), (ii) correspond to Theorems 1.1 and 1.2 of |HM10j . where these statements 
are proved under more general assumptions. In particular, |HM10j does not require the 
uniform upper bound A on the maximal degree. The restriction to the uniformly bounded 
degree sequences implies that the constant R(d n ) used in |HM10j satisfies c < R(d n ) < c _1 
for all n large enough and is therefore immaterial for our purposes. 

Parts (iii), (iv) are taken from Theorems 2.3 and 2.4 of [JL09J. When reading those 
theorems it is useful to realise that the (pj)-distributed random variable D used in |JL09j 
satisfies K[D] = A and that E[D(D — 2)] = XQoo in our notation. 

Remark however that neither [HM10J, or [JL09] consider degree sequences with vertices 
of degree zero, that is with n Q (d n ) > 0. It can however be seen easily, that if n (d n ) does 
not exceed ('n, < 1 (which is implied by the assumptions of the theorem), the vertices 
of degree zero do not have any influence on the existence of the giant cluster, they only 
change the constants A, B in (i), (ii). For (iii), (iv), when no(d n )/n — > po ^ 0, one applies 
the theorem for the modified sequence (d n ) where all vertices of degree are omitted. The 
new degree sequences d n are functions on V n , with n = n(l — po) + o(l). They satisfy 
rii(d n )/n n °°> pi/ (1 — p ) =: pi. Therefore, denoting by the letters with bars the quantities 
related to the distribution (pi), we obtain Q(d n ) = Q(d n ), A = A/(l — p ), f3 = (5/(1 — po), 
g(x) = (g(x) — Po)/(l — Po), and £ = £. This implies that v n = v n and pn = pn, confirming 
that zero-degree vertices have no influence on the asymptotics of the size of the maximal 
connected component. This completes the proof. □ 

4. Random walk estimates 

This section contains estimates on the random walk on random regular multigraphs which 
will be useful later in order to estimate the typical degree sequence of the vacant graph V". 

We start by introducing some notation. We use T d to denote the infinite d-regular tree 
with root 0. For a <i-regular multigraph G = (V, £), a map from T d — > V is said to be a 
covering of G from x £ V, if 4>(0) = x, and for every y G T d , <p maps the d neighbours of y 
m T d to the neighbours of <j)(y) in G, including the multiplicities and the loops. For <i-regular 
multigraphs constructed by the pairing construction this means that the neighbours of y 
are sent by <fi to the vertices which are paired with (<fr(y),i), i G [d]. 

In agreement with our previous notations, denotes the law of the continuous-time 
simple random walk on starting from y. It is important to notice that fixing a covering 
from x, the image by of a random walk in T d with law is distributed as ■ 

For every finite connected A C T d and z G A we now define the escape probabilities 

(4.1) e A (z) = PT[H A = oo], 

which can be calculated explicitly in practical examples using the fact that 

(4-2) Pj d [H = oo} = ^, 

for every neighbour y of 0. This comes from a standard calculation for a one-dimensional 
simple random walk with drift, see for instance |Woe00] . proof of Lemma (1.24). 

We use Zi G T d , i G [d] to denote the neighbours of listed in some predefined order. For 
x G V, let (j) x be a covering of G from x. From now on, if G was obtained by the pairing 
construction, we require that 4> x (zi) is the vertex matched with (x, i). Otherwise (p x can be 
chosen arbitrarily, since our statements will not depend on which particular choice of <p x is 
picked. For any D C [d], we define the sets 

(4.3) B D = {0} U fa : i G D} and B X>D = {x} U {4> x (zi) : % G D}. 
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Figure 1. A covering of a regular graph from x. 



The sets B x n will be used later in the calculation of the degree distribution of the vacant 
set: using an inclusion-exclusion formula, we can express the event {the degree of x in V u 
is k} in terms of events {B xD C V u }, for D C [tfj. 

We first prove a technical lemma describing the graph V after removing the set B x jj. 

Lemma 4.1. For K > 0, we say that the graph G is i^-good, if there is no x G V and 
D C [d] such that: 

(a) B(x,3) is a tree, and 

(b) either V \ B xD is disconnected or diam(V \ B xD ) > Klogn. 
Then, there is K > 1 such that G is K-good, ¥ n ^-a.a.s. 

Proof. To prove this claim we use Lemma 2.14 of [Wor99j which states that P ni( rprobability 
that there exists a (s, j)-separating set in G, that is a set S C V with \S\ = s such that 
G\S contains a component of exactly j vertices, satisfies 

(4.4) f n ,d[G has an (s, j)-separating set] < 3 2+s/d (^-^y (2 ^nHj + s)^ s . 

(This formula is used in |Wor99j to prove that G is P n ^-a.a.s. <i-connected. However, the 
^-connectedness cannot be used directly to prove our claim, since \B(x, 1)| = d + 1.) 

Observe first that it is sufficient to consider D = [d]. Indeed, if x is such that B(x,3) is 
a tree and the graph is connected after removing B x u\, then it is connected after removing 
B x D for any D C [d]. In addition, obviously diam(G \ B x D ) < diam(G \ B x ^) + 2. 

Let now D — [d], that is B x D = B(x, 1). Assume that G\B(x, 1) is disconnected. Then 
removing the set B(x, 1) \ {x} from G, divides the graph into at least three components. 
One of them is {x}, and the other ones are are contained in G \ B(x, 1). Since we require 
that B(x, 3) is a tree and since G is <i-regular, the size of these other components is at least 



4. We thus apply (4.4) with s = d and j > 4: The P n]( rprobability that there is x G V such 
that B(x, 3) is a tree and G \ B(x, 1) is disconnected is bounded from above by 

n/2 

(4.5) ^^P ni d[G has an (s, j)-separating set]. 

i=4 

The largest term in this sum, corresponding to j = 4, is 0{n 4 ~^ d ) = o(l). All the remaining 
terms are much smaller. Actually, as in the proof of Theorem 2.10 of |Wor99j . it can be 



shown that (4.5) tends to as n — » oo. Hence, P n ^-a.a.s. there is no x G V such that B(x, 3) 
is a tree and G \ B(x, 1) is disconnected. 

To bound the diameter of G\B(0, 1) we use the fact that P^-a.a.s. diam(G) < 2 logj^ n, 
sec [BF82J. We may also assume that G \ B(x, 1) is connected, which as we proved occurs 



We now claim that 

„s removing one vertex v of degree d in an arbitrary graph while keeping it 
connected can increase the diameter of the graph at most by a factor of 3 d-1 . 

To prove this claim we first consider the removal of an edge: Removing one edge e from a 
graph while keeping in connected can increase the diameter of the graph at most by factor 3. 
To see this it is sufficient to consider the shortest path \i in G \ {e} connecting the vertices 
of e (such path must exist since G\ {e} is connected, and cannot be longer than 2diamG), 
and to replace the edge e by the path \x in every geodesic of G that contains e. 

Having understood the removal of edges, we can analyse the removal of a vertex v. We 
first remove all but one of the edges of G incident to v, in this procedure the diameter of 
the graph is multiplied by at most 3 d_1 . Removing v together with the last edge linking it 
to G \ {v } yields the claim (4.6). 

The claim (4.6) and \B(x, 1)| < d + 1 then imply that P ni rf-a.a.s 

(4.7) 



diam(G \ B(x, 1)) < 3 { ^ 1)(d+1) diamG < 3 {d - 1)(d+1) • 21og d _ 1 
This completes the proof of the lemma. 



n. 



□ 



We now start controlling how the random walk visits the sets B X £, defined in (4.3). In 
the lemma below, we show that the probability of escaping from B x d for a large time can 
be approximated by the escape probability on the infinite tree, defined in (4.1). 

Lemma 4.2. For every x G V n , D C [d] and i G D 



(4.8) 

and 

(4.9) 



E, 



n,d 



P n AHb x , d <log 2 n] < 



< 



clog n 



n 



clog n 



n 



Proof. To simplify the notation we write B, B, z and <j)(z) for Bp, M Xt u, Zi and <f> x (zi 
the fact that (j) x maps a random walk on to a random walk on G, we can write 



Using 



(4.10) 



Z) [H B < log 2 n]=Pj [H^ [B) <\og 2 n] 



< loe 2 n, Hu > loe 2 n] 



Therefore, the left-hand side of (4.8|) can be bounded from above by 
PT [H M >\og 2 n]- 



(4.11) 



+ E. 



n.d 



Using the Markov property at time log n and the definition of 



< log n] 

n(z), the first term equals 



PT [log 2 n < H m < oo] 



(4.12) 



< P 



d(0,X Y 2 n ) < 



d-2 
~2d 



log n 



+ sup Pi 

u:d(u,0)> l°g 2 n 



H < oo 



Since d(X t , 0) under Pj d is a random walk on N with expected drift given by (d — 2)/d, 
both terms above are bounded by cexp{— d log 2 n}. 

To bound the second term in (4.11 ), note that if at time t the random walk on started 
from z visits a point in (f)~ l (B X £>) \ B, then the trajectory of the image walk on G together 
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with B(<p(z),2) contains a cycle in G. Therefore, denoting by G\a the subgraph of G 
generated by A C V n , 



(4.13) 



E 



n,d 



< log 2 n] 



< P 



n,d 



G\ 



B{X Q ,2)U{X a :s<log 2 n] 



contains a cycle 



Taking care of the possibility that the continuous-time rando m w alk makes more than 
21og 2 n steps before time log 2 n, using the notation from Section 2^3 this is bounded from 
above by 



(4.14) 



PT[N lo ^ n >2\og 2 n}+P n jG\ B{X(h2 



)U{X i: j<21og 2 n} 



n i contains a cycle 



The random variable iV log 2 n has Poisson distribution with mean log n, therefore the first 

term in (4.14) is smaller than ce~ c ' log2 ™. To bound the second term, observe that the con- 
sidered subgraph can be constructed inductively by the following variant of the construction 
from Section l2~2l 

(1) Choose Xq uniformly at random in V n . Use the pairing construction of Section 2.2 
to construct the set B(X , 2). This requires creating at most d + d(d — 1) pairs. 

(2) Let X\ be a uniformly chosen neighbour of Xq. (All neighbours of Xq are known 
from the first step.) 

(3) Repeat for all i G {2, . . . , 2 log 2 n} the following steps: 

(a) Choose uniformly in [of], independently of the previous randomness. 

(b) If the half-edge (X^i, Zj) is not yet matched, then match it with an half-edge 
chosen uniformly among the remaining half-edges, as in the pairing construction. 

(c) Let Xi be the vertex that is matched with the half-edge (Xj_ 1; Z{). 

The probability that a cycle is created in the step (1) is easily bounded by c/n. It is 
not possible to create any cycle in the step (2). In the step (3) the cycle is created only 
when the half-edge (Xj_ 1; Zj) is not yet matched (otherwise we do not add any new edge 
to the subgraph) and when the vertex Xi was already 'visited by the algorithm', that is it 
has some matched half-edges from the previous steps. Since the algorithm visits at most 
d + d(d — 1) + 21og 2 n vertices, the probability to match (Xj_i, Z{) with an already visited 
vertex is smaller than cn _1 log 2 n. Therefore the probability to create the cycle in the 
step (3) is at most cn _1 log 4 n. In consequence, the second term in (4.11) is smaller than 
c'log n _|_ cn - l \ g A n which establishes (4.8). 



ce 



In order to prove (4.9), observe that 
(4.15) P n A H B < log 2 n] < P M [iV log2n > 21og 2 n] + £ P n „ 



2 log 2 n 



Xi G B 



i=0 



Since under P n> d, the random vertex Xi is uniformly distributed in V n , P n< d[Xi G B\ < 



n 1 \B\ <{d+ l)n 1 . Hence, 
(4.16) P n A H B<\og 2 n}<c'e 



clog n 



+ 2(d + l)^ 1 log 2 n < en' 1 log 2 



n. 



This completes the proof of Lemma 4.2 



□ 



The previous lemma has a simple corollary. 
Corollary 4.3. Fix D C [d] and i G D, then both F n ^-a.a.s. and F n ^-a.a.s. 



(4.17) 



{xeV n : \Pl {zi) [H Bx>D >\og' 



n\ 
ll 



> n" 1/2 } < (log 5 n)n 



1/2 



and 
(4.18) 



{x e V n : P G [H Bx D < \og 2 n] > n' 1 ' 2 } < (log 3 ™)™ 1 / 2 . 



Proof. Note that the complements of the above events are respectively contained in the 



events ^ 



\pG 



< lOE 



&)l > log a n and E x ev n P l H B x , D < logfn] > 



log n. Thus, the P nj d-a.a.s. statements of the corollary are implied by the Markov inequality 
and Lemma 4.2 The P n(t j-a.a.s. statements then follow using the remark below (2.3). □ 

In what follows we will make use of the quasi- stationary distribution with respect to a set 
B C V. The usual definition of this distribution is given in the Appendix, see (A.l). The 
quasi-stationary distribution can be thought as the asymptotic distribution of the position 
of a random walk conditioned not to visit B. To proceed, we will need the following lemma 
which gives the rate of convergence of such conditioned walk towards the quasi-stationary 
distribution. 

Lemma 4.4. Let t n = log 2 n. Then, ¥ n t i-a.a.s., for any x such that B(x,3) is a tree and 
for any connected A C B(x, 1) 



(4.19) 



sup 

x,yeV\A 



< ce 



-c' log 2 n 



P^[Y tn =y\H A >t n }-a A (y) 

x,y&V\A 

Proof. By Lemma 4.1, V\A is P^-a.a.s. connected and has diameter smaller than K\ogn. 
To prove the claim of the lemma we are going to make use of Lemmas A.l and A. 2 , which 
are general results on Markov chains presented in the Appendix. We use the notation 
introduced there, in particular < \± < \^ < ■ ■ ■ stand for the eigenvalues of the (— C A ) 
where C A is the generator of the random walk killed on hitting A. 

From (2.10), Ac > a > 0, P„ 5 d-a.a.s. Therefore, using Lemma A.l, we obtain that 

(4.20) X A - Xf > a - ' 



E[Hj 



Pn.d a.a.s. 



To obtain an upper bound on E[Ha] \ we use Proposition 3.1 and (2.10) of [CTW11 
Using this proposition with A = A and C = V \ A, we obtain 

(4-21) < cir(A) < C 



This implies that P n ^-a.a.s. 



E[H A ] 
\ A >c>0. 



n 



Using the above fact together with Lemma A. 2 (observing that tt(x) = 
graphs) we obtain that P^-a.a.s. 

f 3 / 2 exp{— clog 2 n} 



1/n on regular 



(4.22 



sup 

x,y£V\A 



P^[Y tn =y\H A >t n }-a A (y) < 



cn 



inf ze v\A cta(z) 

We have to bound the infimum in the denominator. For this, take z 6 V \ A and any 



t > 0. By reversibility, for any z' E V \ A, 



(4.23) 



P$[X t = z\H A >t}= Pf[X t = z'\H A > t] 



pj?[H A >ty 



In order to bound the above ratio, note that 



(4.24) 



P?[H A >t]> P?{H Z , < H A , H A o 9 H , >t] = P^[H Z , < H A }P^[H A > t}. 



G\ 



As P n> d-a.a.s. the graph induced by V \ A has diameter at most K log n, we can find a 
path 7 with length at most K log n, connecting z and z' and not passing through A. This 
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gives us that P Z [H Z > < H A ] > d~ clogn > cn~ c ' . From Lemma A_2 lim^oo P w [X t = z\ H A > 
t] = &a{z) uniformly for all w, z G V \ A. Therefore, taking the limit t — > oo in (4.23), 
&a(z) > c(TA{z')n" c . Together with the fact that a a is a probability measure this yields 



(4.25) 



inf o 'a(z) > cn 

z£V\A 



Using the above result with (4.22) finishes the proof of Lemma 4.4 



□ 



5. Degree sequence of the vacant graph 

We are now in position to estimate the typical degree sequence of the vacant graph Vq 
under P nj rf. At this point it is instructive to mention the relation between this set and the 
random interlacements process on the ci-regular tree T d . 

The model of random interlacements on transient weighted graphs is constructed in 
|Tei09j . see also |SznlOj for the original construction of the model in the particular case 
of Z d , d > 3. Random interlacement on T d can be understood as a measure Q u on the space 
{0, 1} T which samples a random subset Vj d of T d (called the vacant set left by random 
interlacements at level u) characterised by the following: 

(5.1) Q U [K C V£] = exp{-u cap(iC)}, for every finite K C T d , 



where cap(K) = ^2 xEK e^-(x) for ex as in (4.1), cf. (2.27) of [Tei09j. 



Intuitively speaking, the vacant set of the random interlacement Vj d gives the asymptotic 



local picture of Vq under ~P n ^ as n tends to infinity, see Proposition 6.3 of [CTWllj . An 
important fact about random interlacements on T d is that the law of the vacant cluster of the 
root G T d under Q u is the same as the law of a certain (inhomogeneous) Galton- Watson 
tree. More precisely, the probability that is vacant equals e - nca p( ) = eX p|_ M ^| j. Given 
that G Vj d , the offspring distribution of the root is binomial with parameters d and 

(5.2) p u = exp <^ - u- 



d(d-l) 

and the offspring distribution of all remaining individuals is binomial with parameters d — 1 
and p u . Using this characterisation it is easy to compute the probability that the degree of 
the root in Vj d is i, i = 0, . . . , d. 

:= Q u [degree of in Vj d equals i] 

Q u [0 is vacant and has exactly i offsprings] 



(5.3) 

d 



a— 1 



We now explain the relation between df and the degree sequence of V G . This relation 
was a lready ob tained in a weaker form in Theorem 3 of |CF10] and could also be extrac ted 
from jCTWllj Proposition 6.3. However the finer control of errors obtained in Theorem 
is crucial if one wants to stay inside the critical window. 



5.1 



Recall from Section [3TT] that T> u denotes the degree sequence of the vacant graph V" and 
that, for any degree sequence d, rii(d) denotes the number of vertices with degree i in d. 

Theorem 5.1. For every u > and every i G {0, . . . , d}, 

(5.4) \E G [m(V u )} -n<%\ < c(\og 5 n)n 1/2 , f n4 -a.a.s. 
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Proof. We fix x G V n , D C [d], and recall from Section [4] the definitions of the covering map 
0, of B xD , M D and of Zj. To simplify the notation we use a and B as shorthand for o BxD 
and B X) d. We first estimate the probability that B C V. Using the Markov property and 



Lemma 4.4, we obtain that P^-a.a.s 
P G [H B > un] - exp j - 



tin 



(5.5) 



< 



E G [H B 

P G \H B > t n ]E G [Pg n [H B >un- t n ]\H B > t n ] - exp j 
P G [H B > t n ]P G [H B >un- t n ] - exp { 



un 



un 



E?[H 



ce 



E g [Hb] 

— c' log 2 n 



Under the measure P a , the random variable H B is exponentially distributed, see for instance 
[AB93] below (12). Hence, using that e~ l is 1-Lipschitz for t > 0, 



(5.6) 



P G [H B >un-t, 



exp 



un 



E a [Hi 



< 



E G \H B \ 



Therefore, ( 5.5[ ) becomes 



(5.7) 



P G [H B > un] - exp j 



un 



E?[Hb 



Let K 00( j C V be the set of vertices x G V satisfying 



(5i 



B(x, 2) is a tree, and for every D C [d] and i £ D, 
\P^)\Hb x , d > tn] - e B (0(^))| < n-V2 and P G [H BxD < t n ] < n^\ 



By Corollary 4.3 above, and by Remark 1.4 and Lemma 6.1 of jCTWllj . the complement 

\V \ U good | < cn^Qogn) 5 , P M -a.a.s. 



of V^ 00 d is very small, 



(5.9) 



By Lemma 2 of |AB93] and p~2~Ij) 



(5.10) 



£7?[if B ] 1 < £ G [if B ] 1 < cn~\ 
Therefore for x G V^d, (5.7) becomes 



(5.11) 



P G [Hs>un]-^{-^-\ 



< en- 1 ' 2 . 



Our next step is to obtain an estimate on E a [H B ] for x G V^ ood . To this aim we 'collapse' 
the set B X} d into one point b, and define a new Markov chain whose distribution is denoted 
P and which is characterised by its transition rates p xy , 



(5.12) 




= 2W, i£w,w'^b, 



It is easy to check that 7r = + J2 x ^b^^) * s a reversible distribution for this chain. 

Therefore, 



(5.13) 



n 

W\ 



Eb[Hb] = E b [H b , H b < t n ] + P b [Hb > t n ]Eb [Ex tn [Hb — t n ]\H h > t n ] . 
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By ( |5l2) , P 6 [if 6 > * n ] = (fip 1 J2 y eB P y[ H B > *tj- Therefore, using Lemma |4.4| and ( [5~l3] ), 



(5.14) 



£P y G [F B >t n ]i^] 
yes 



n 



< c£„ + cexp{— c't n } < ct Ti 



Using (5.10), this yields the following estimate on E G [H B ] 
(5.15) ■ » ^ <*.■ 



< 



clog n 



n 



We are now in position to give our final estimate on P g [Hb > un]. By the triangle 
inequality, for x G V goo d, 



P G [H B > un] - exp | - n^e B (y)| 

j/eB 

P G [F B >H-exp{-^-^ 



(5.16) 



+ 



+ 



exp { - iim } - exp { - n S p " ^ > tn] } I 



exp | - u P y\ H B x , D > t n ] } - exp | - u ^ e M } 

yeB x , D yen 



< cri 



-1/2 



where for the last inequality we used the estimates (5.11), (5.15) and (5.8). 

If x G V^ 00 d, all its neighbours are distinct. We can then use the inclusion-exclusion 
formula to write 



(5.17) 



P G [V u (x)=i]= P G [V u r\B(x,l) = B x , c ] 

Cc[d],\C\=i 

£ (-i) |DHC| p c [ff B „ >H' 



E 

C<Z[d},\C\=i Dd[d],C<ZD 



Using (|5.16|), we obtain that 
(5.18 



P G [V^x) = i ]- £ £ (-ir\-Mex V {-uJ2e®M} 

Cc\d],\C\=i Dc\d],CcD y& D 



< cri 



-1/2 



From (4.1), (4. 2), it is not difficult to deduce that for y G B Xy o \ {x}, 



(5.19) 



d- 1 d-2 , . , d- \D\ d-2 

and e Bn m = 



d d — 1 



Inserting this into (5.18) leads to 
(5.20) 



P G [V u {x) = i] 



3=1 



d — i 



exp < — u 



d d — 1 



d-2/. d-l d-j 



d-l\ 3 ' d 



+ 



d 



< en' 1 / 2 . 



A simple computation implies that the leading term in the last formula equals df (see (5.3)). 
Therefore, P^^-a.a.s., uniformly for x G V goo d, 

(5.21) 



\P G [V U ' 



x) — i] — d"| < cn 
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-1/2 



The claim of Theorem 5.1 then follows by summing this relation over x G V^ 00 d and using 



We now prove the concentration of rii(T> u ) around its mean. 

Theorem 5.2. Let G be a d-regular (multi) graph on n vertices satisfying \q > a > 0. 
Then, for every e G (0, |), and for every i G {0, . . . , d}, 

(5.22) P G [\ ni {V u ) - E G [ni{V u )} \ > n l/2+£ ] < c £ e~ cne . 

To prove this lemma we use the following concentration theorem for (not necessarily 
independent) random variables that we learnt from |McD98j . Consider a sequence W = 
(W±, . . . , Wm) of random variables, all taking values in some space A. Let / : A M — > R be 
a bounded function. For k G {1, . . . , M} and yx, . . . , yu-i G A^ 1 we define 

n(yi, ■ ■ - ,2/fc-i) 

(5-23) = gup \E[f(W)\W k = y,W t = y t Vi<k}-E[f(W)\W k = y',W t = y t Vi<k}\ 

y,y'eA 

and set 

M 

(5.24) R 2 = sup | 5^rjj(yi, . . . ,y*_i) : yi, .. . ,2/m-i G a\. 

k=l 

Lemma 5.3 (Theorem 3.7 of |McD98| ). Let W = (W u . . . , W M ) be as above. Then 

(5.25) W(W) - Ef(W)\ >t]< 2e~ 2t2/R2 . 



Proof of Th eorem 5.2, To apply Lemma 5.3, we need the following co nstr uction similar to 
Section 4 of jCTWllj . Let £ = n £ for e from the statement of Theorem 5.2 On an auxiliary 



probability space (Q,Q), define (Z % ,i G N) to be a collection of i.i.d. uniformly chosen 
vertices of G. Given the collection (Zi), let (Yi : i > 1) be (conditionally) independent 
family of elements of D(\0, £], G ) such that Y % is distributed according to the random walk 



bridge P z Li z % (see Section 2.3 for the definition). We define X G D([0, oo), G) to be the 



z i ~ 1 ,z i 
concatenation of Y l, s, 

(5.26) X(t) = Y\t- (i- 1)£), when (i - 1)£ < t < ii. 

We use V G to denote the distribution of X on D([0, oo), G), V G = Q o X" 1 . V G > T stands 
for its restriction to D([0,T],G). The measure r p G -> un approximates well F G,m if £ is large 
enough as follows from the next lemma whose proof is postponed to the end of this section. 

Lemma 5.4. Assume that Xq > a and £ = n £ . Then there exist constant c a ^ £ and c' a£ such 
that for every u > and all n satisfying ne~ £a < 1/2, p G ' un and P G > un are equivalent and 

(5.27) 



dP G > 



1 



< c'ue c ^ 



To be able to apply Lemma [5.3[ more precisely to estimate the functions r k , we do not 
want | RanF'l to be too large. Therefore, we define Y % C V to be the set of first 2£ vertices 
visited by Y'\ 

(5.28) Y i = {Y?,t<(T 2t (y i )Ae)}, 

where Tk(Y l ) denotes the time of fc-th step Y l (defined to be infinite if Y l makes less than 
k steps). Obviously Y % C RanF*. On the other hand, it can be proved as in Lemma 4.2 of 
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CTWllj . that 
(5.29) 



Q(RanY i ^ Y l ) < sup P% £ [N £ > 2£\ < ce~ c ' 1 . 



x,y£V 



(Remark that Y % has the law of the random walk bridge, and thus (5.29) is not just a direct 
consequence of large deviation estimate for a Poisson random variable.) 

We may now prove Theorem 5.2 Set m = \un/£\ and u' = m£/n. Let N be the number 
of steps of X between u'n and un. Since un — u'n < £, by properties of Poisson random 



variables, P G [N > 

(5.30) 

and 



< e 



-cl 



Between, u'n and un the walk visits at most N sites, therefore 
\ ni (V u )-n t (V u ')\ < (d+l)N, 



(5.31) \E G [ ni (V u )} - E [ni(V u )]| < (d + 1)E G [N] < (d+l)t 
Thus, for n > c a>e , 

P G [\m{V u ) - E G [ ni {V u )\\ > n l ' 2+£ ] 

(5.32) <P G [\n l iV u ')-E G [n i iV u ')} \ > n 1/2+£ - 3(d + 1)£] +P G [N>2£]. 
< P G [\ni{V u ') - E G [n i (D u ')}\ > \n l / 2+£ ] + e 



-cl 



Using Lemma |5.4[ denoting by E G the expectation corresponding to V 



G 



(5.34) 



(5.33) \E G [ni{V u )] - E G \n { {V u )]| < c a , € nu'e- c ^ e , 

and thus, for n > c Qje , 

P G [\n l {V u ')-E G [n l {V u ')}\ > \n l ' 2+e ] 

< V G [\n i (V u ')-S G [n t (V u ')} \ > {n 1 ' 2 ^] + c a>t u' e'^ . 

Observe that under P G , V"' = V \ Uj< m RanF l . Let V be the vacant set left by F l 's, 
V = V \ Uj< m F l , and denote by T> the degree sequence of the graph with set of vertices V 
and edge set {{x,y} G £ : x,y G V}, cf. (3.1). By (5.29), we have then 

(5.35) Q[V U ' ^V]< cme- c ' £ < ce" c ' 7 . 

Therefore, for n > c„, 



(5.36) 
and thus 

(5.37) 



£ G [m(V u ')] - Q[rn(p)]\ < nQ[V u ' ^ V] < ce 



-c'£ 



V G [\n l (V u ')-S G [n t (V u ')}\>y/ 2+£ ] 

<Q[\ni{V) - Q[ni(V)]\ > In 1 ^] + ce~ c ' £ . 



We now apply Lemma 5.3 with M = m, W, t = Y l , f = ni(T>), and A being the set of 

(yi, 



subsets of V with at most 2£ elements. Writing y k = (y 1; 
and Yk = (Yi, . . . , Yk), we claim that 

(5.38) r k {y k _ x ) = sup \Q[ ni {V)\Y k = y k ] - Q[n t (V)\Y k = y' k }\ < 2(d+l)£. 

y,y'eA 

Indeed, by conditioning also on the values of Y k+2 , . . . , Y m , we observe that the difference 

(5.39) \Q[n t {V)\Y k = y k , Y k+2 , . . . , Y m ] - Q[n t (V) \Y k = y' k , Y k+2 , . . . , Y m ] \ 



cannot be larger than (d+l)(\Y k \ + \Y k+1 \) < 2(d+l)£. The inequality ( |5.38D then follows 
by integrating over Y k+2 , . . . , Y m . 
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From (5.38) it follows that we can apply Lemma 5.3 with R 2 = m{d + 1) 2 £ 2 
yielding 

Q[\ni(V) - Q[ni(V)]\ > \n 1/2+£ ] < ce~ cn w ^ nl+£ 



a i 



l+E 



(5.40) 



< ce" 



This, together with (5.32), (5.34) and (5.37) completes the proof of Theorem 5.2 



□ 



Proof of Lemma 5. J/.. Let v! be the smallest number greater or equal to u, such that u'n is 
an integer multiple of £, and set m = u'n/£. Let further A be an arbitrary J-" uri -measurable 
subset of -D([0, u'n], V). Since P G ' un and , p G ' un are the restrictions of the measures 

pG,u'n 

and "P G ' un to D([0,un],V), it is sufficient to prove the lemma with u replaced by u' . To 
this end we write 



(5.4i; 



P 



G,u'n 



[A] 



P G > u ' n [A\X ie = x u 0<i< m]P G > u ' n [X ie = Xi , < i < m]. 



X0,...,Xm£V 

By the Markov property 



m— 1 



(5.42) 



P G ' u ' n [X t£ = Xt ,0<t<m] = vr(x ) n P* k [X e = x k+1 ]. 



k=0 



The construction of the measure V G,u ' n implies that 

V G ' u ' n [A\X i£ = Xi ,0<i<m] = P G ' u ' n [A\X ie = Xi,0<i<m] 

(5.43) 



V G ' un [X, 



it 



Xi, < i < 



2m] = J] 



7T(X k . 



k=0 



Comparing (5.42) and (5.43), it remains to control the ratio -Pf[X^ = y]/ir(y). However, by 
(2.9) and the assumption of the lemma, \P G ' l [Xi = y]/Tr(y) — 1| < ne~ a£ . This leads to 



-pG,u'n r Al 

(5.44) (1 - ne a T < pG^] ^ ^ + ne ~^ m 

Since I = n e and ne~ £a < | by the assumptions of the lemma, it immediately follows 
that P G ' u ' n and V G ' u ' n are equivalent. A change of constants accommodating the terms 
polynomial in n then completes the proof. □ 



6. Proofs of Theorems 11.11 and 11.21 

We now have all tools that we need to show all the main results of this paper. As a direct 



consequence of Theorems 5.1, 5.2 and the fact (2.10), we get that P nj rf-a.a.s. 

\ni(V u ) -nd%\< cn 1/2 log 5 n, for all < i < d, 



(6.1) 



where is defined in (5.3). Hence, Pji,d~&-&-s, lim n _j. 00 n ^77^(2}") — d^ . The constant 



Q{V U ) (see (3.6)) can be written as 
(6.2) 



Thus, for n > c u , on the event in ( |6.1 ) we have 



(6.3) 



Q(V U ) - ( Pu (d - 1) - 1) 
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< c„n 1 ' 2 lo£ 



n. 



The value given by (1.2) satisfies p Uit (d — 1) — 1 = 0. Therefore, when u n — > u+, we 
obtain by expanding the exponential in the definition (5.2) of p u around u±, 



(6.4) 



Q(V U ") - (u+ - Ur. 



(d-2) s 



d{d -I) 



< c({u+-u n ) 2 + n~ 1/2 log 5 n). 



This implies that when u n is in the critical window of The orem 1.1 , that is \n l ^{u^—u n ) \ < 
A, then Q(V Un ) is in the critical window of Theorem pjjl that is n l ^\Q{V Un )\ < A', P n>d - 



a.a.s. Theorem [T7T] then follows directly from Theorem 3.2 i) together with Proposition 3.1 
and the remark following (2.3). 



Very similar reasoning apply when proving Theorems 1.2 and 1.3 We should only identify 



the constants of Theorem |3.2| Easy computations give 
(6.5) 

and thus 



X = e - u *=id Pu , (3 = e~ u ^d(d- \){d-2)pl, 



(6.6) 



2n\ 2 f3- 1 Q(V Un ) = 2n(u+ - u r 



d-2 
Jd^l) 



■e-"*^i(l + o(l)). 



Replacing u± — u n by u n n 1 / 3 , Theorem 1.2 a) follows. It can also be seen that >Jn/Q {V> Ur - 
is of order n 2 ^ 3 u n 1 ^ 2 , implying Theorem 1.2 b). 



Finally to identify p of Theorem 1.3 We observe that g(x) of Theorem 3.2 is given by 



(6.7) 



g(x) 



J2d^ = e- u ^(x Pu + (l- Pu )Y 



i=0 



X, 



After few simplifications, £ of Theorem 3.2 is the unique solution in (0, 1) of the equation 

(6.8) (xp u + (l-pu)) d - 1 
and p is given by 

(6.9) p=l-g(£). 



This completes the proofs of all three main theorems. 



Remark 6.1. (1) Theorem 5.1 raises the question of what is the right magnitude of devi- 
ations in E G {rii(T) u )] under the law P nj d? If indeed it is of order n 1 ^ 2 (without power-log 
corrections), then it would be interesting to investigate whether this quantity satisfies a 
central limit theorem when properly rescaled. 

and Theorems 



(2) As established in Proposition 3.1 



5.1 and 5.2, we can reduce the study 



of V u to questions on the behaviour of random graphs with prescribed degree sequences. 
Although the results in |MR95l IJL091 IHMlOj provide very fine information about such 
graphs, there are several questions concerning them which are still open. For instance, one 
could give a better description of the geometry of their critical components, their diameters, 
spectral gaps, etc. 

(3) It is interesting to notice that the statements (1.3) and (1.4) were established in 



CTWll] for the vacant set left by random walk on other sequences of graphs, such as large 
girth expanders. Is it possible to extend the results of the current paper on the critical 
behaviour of V u to this more general setting? 
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Appendix A. Properties of the quasi-stationary distribution 

We establish here few results for arbitrary reversible irreducible continuous-time Markov 
chains on a finite state space. These results are natural but we have not found any suitable 
formulation in the literature. 

Let V be a finite set and let £ be the generator of a reversible irreducible continuous- 
time Markov chain X on V, and let ir(x) be its invariant measure. We use (/, g) to denote 
the usual scalar product on L 2 (V, ir), (f,g) = J2 xe v fi^d^) 71 ^)- The operator — C is 
symmetric in L 2 (V, ir) and has real eigenvalues = Ai < A2 < • • • < Aiyi and corresponding 
orthonormal eigenvectors Vi, . . . , v\v\- 

For B C V, we use L B to denote the generator of the Markov chain X killed on hitting 

B, viewed as an operator on L 2 (V \ B,tt\v\b)- Let < Af < Af < ^\v\b\i an d v i > 

. . . , t)im B i denote the eigenvalues and eigenvectors of — C B . 

The quasi-stationary distribution a B is related to the eigenvector of — C B corresponding 
to Af and is given by 

{AA) aB{y) - l^ir " Wv' 

with 1 denoting the constant one function. Inverting this relation we get 

(a.2) »fw=^( e ^r /2 . 

1 ir(x) V ^ tt(x) J 

V ' xdV\B v ' 

Lemma A.l. For every B C V , 

(A.3) A? _ Af > A2 __l_. 

Proof of Lemma A.l\ Since C B can be viewed as a sub-matrix of C, by the eigenvalue in- 
terlacing inequality (cf. |Hae95] . Corollary 2.2), we have X B > A2. On the other hand, by 
|AB93j Lemma 2 and the paragraph following equation (12), 

(A ' 4) A ? = OT"W 

Combining these two inequalities we obtain Lemma [A.l □ 



Lemma A.2. Suppose that for t > and e G (0, 1/2) 

°b{x) \ 2 o- B {x) 



(A.5) e-«^-W\V\B\( sup ^L) < e inf 

\xev\B a/7t(x)^ xey\B a/7t(x) 

Then, 

(A.6) sup |P x [X t = y\Hs >t\- a B {y)\ < As. 

x,y£V\B 

Proof. In the proof we will only use the eigenvalues and eigenvectors of —C B , therefore 
we omit the superscript B from the notation. Similarly, we write a for <7b and define 
m = \V\B\. By the usual spectral decomposition formula, 



P x [X t = y,H B >t) = (e tcB 5 y )(x) = J2 e ~ Xkt Mx)(vk,5 y ), 

(A.7) 

P X [H B >t] = (e tcB l)(x) = J2 e ~ Xkt Mx)(v k ,l), 



fc=i 
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where 5 y is the indicator function of y. For / e L 2 (n), define ipf = Y^k=2 e ^ i v k> f) v k- 
Then e tcB f = e~ Xlt (vi(vi, f) +i/)f), and by Pythagoras' theorem 

(A.8) HV/lbw < e- (A - Al) 1/|| L2W . 

Using this notation and the definition of the conditional probability, 

vi{x){vi,5 y ) +4>s y (x) 



(A.9) 



P x [X t = y\H B >t] 



Applying (A.l) we get after an easy algebra 
(A.10) P x [X t = y\H B >t]-o(y) 



vi(x){vi, 1) + ipi(x) 



Let / stand either for 5 y or 1. Then H/Hz^) < 1, which directly implies \ipf(x)\ < 



7r(x) 1 ^ 2 e ( A2 Al * ) *. From (A.8), (A. 2), using the assumption (A. 5), we obtain 



(A.ll) 



ip f {x) 



-(A 3 -Ai)t 



< 



"!-■'-' e-( A2 ~ Al )'m sup z ( -^=) 2 



< 



inf. 



\Ao) 



< e. 



Using e < 1/2, this implies that the absolute value of (A. 10) can be bounded from above 
by 

2(|^(a;)i + iVi(x)i) 



(A.12) 



< 4e. 



This completes the proof of Lemma A. 2 



□ 
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